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Abstract 

We explore the stability of the variance and skewness of the cosmic gravitational convergence 
field, using two different approaches: first we simulate a whole MEGACAM survey (100 
sq. degrees). The reconstructed mass map, obtained from a shear map recovered by the 
usual method via a maximum likelihood method, shows that the actual state-of-the-art data 
analysis methods can accurately measure weak-lensing statistics at angular scales ranging 
from 2.5' to 25'. We looked also at the influence of a varying signal-to- noise ratio over the 
shear map (due to local variations of source density) on the mass reconstruction, by means 
of Monte-Carlo simulation. The detectable effect at small scales can easily be corrected-for 
in most of the relevant cases. These results enhance the confidence in the capability of future 
large surveys to measure accurately cosmologically interesting quantities. 

1 Introduction 

The cosmic shear effect has recently been detected (see Bacon et al. 2000, Wittman et 
al. 2000, Kaiser et al. 2000, van Waerbeke et al. 2000 and these proceedings) and the 
shear variance has been measured in a range of scales from 30" to 10'. However, theoretical 
investigation of lensing statistics (Bernardeau et al 1997) show that the measurement of the 
third order moment of the convergence k is necessary to break the degeneracy between Q,m 
and cTg (for symmetry reasons, the third order moment of the shear is null) . This implies to 
proceed to a full reconstruction of the convergence map. 

Knowing that the skewness of the convergence field is a statistic much more difficult to 
measure than the variance, this raises a few questions: (1) do the measurement errors preserve 
the correlation between cells? (2) How the measurement errors of the shear propagate into 
the convergence reconstruction? (3) How to correctly estimate the noise due to the intrinsic 
ellipticity of galaxies in an observational context? (4) How clustering of the galaxies (leading 
to a spatially varying S/N ratio) affects this reconstruction? 

The points 1) and 2) have been adressed by van Waerbeke et al. (1999), considering only 
the noise coming from intrinsic ellipticity. We tried a different approach: we performed two 
kinds of simulations to partly address all those problems. 

2 The simulated lensing survey 

We performed the simulation of a large (100 deg^) compact survey whose characteristics are 
similar to those of the future MEGACAM dedicated survey at CFHT. The simulation is 
generated in 6 steps: 

First we have created the synthetic catalog of the whole survey (in an open CDM Universe) , 
using the SkyStuff " program. We then modified the relevant galaxy parameters (axis ratio, 
luminosity, etc.) accordingly to 2-dimensional simulations of dark matter with an angular 
resolution of 2.5' (see van Waerbeke et al (1999) for details). We generated in this manner 
numerous sub-catalogs (40x40) and used the SkyMaker^ software to create the synthetic 
images of the galaxies lensed by this 2D simulation. These images have been analyzed with 
/MC^r[[ and we used SExtractor to clean the catalogs of spurious detections (and to select 



"Sky Stuff and SkyMaker are freely available a t 
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Size 


100 sq. degrees 


Ca-talog initicil density 


~ 31 galaxies/arcmin"^ 


Density of detected objects 


~ 22 galaxies/arcmin*^ 


CcLta,log fecitures 


I'igals t [lijZ'^.OJ, IjiVivy exL. law 


Image pixel size 


0.2" 


Seeing FWHM 


0.7" 


PSF anisotropy 


0.05" rms/axis (tracking) 


Sky SB, read-out noise 


20 mag/arcsec^ (I Band), 5e"/ADU 


Object selection 


FLAG = 0, jei.al < 1 


Stars 


m, G [10,21.5], 60 usables/image 



Table 1: Characteristics of the simulated survey: the catalog is obtained from a Schechter function and a PDE prescription, and designed 
to reproduce observations. The images are simulation of CCD frames taken at CFHT. 



the objets with uncorrupted photometry as indicated by the FLAG keyword). The measured 
shear was then used to reconstruct the k map. 

The characteristics of the survey are listed in table |^. As we have only 2-dimensional dark 
matter simulations for a lOx 10 degrees map, all sources in a given area are affected the same 
way, regardless of their redshift. Considering the wide variety of shear values and galaxy 
parameters, this should not weaken our conclusions. The shear is estimated by the quantities 
where the P'^'s are components of the shear susceptibility tensor and <> 



<ei> 



and 



<e2> 



means a spatial average. 

To correct for the noise coming from intrinsic ellipticity, we have simply calculated 
Since the theoretical noise is defined from the complex ellipticity, 
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which is not equal to the KSB-defined ellipticity (which is itself filtered by Gaussians win- 
dows, etc.) it is interesting to get such a simple estimator. Nevertheless, this estimation is in 
very good agreement with the noise one can measure directly on the power spectrum of I7I, 
indicating that we obtain a white noise. We applied the results of Bernardeau et al (1997) 
to correct for the cosmic variance, using a 'sample size' equal to 6 times our largest scale of 
measure, i.e. 2.5 degree. The variance and skewness of k that we obtain as a function of the 
filtering scale are shown in fig 1. 

Although the simulated images still have a quality superior of those of real data (e.g. no 
strong PSF anisotropy - we still performed the anisotropy correction -, no big diffraction 
spikes, no bad columns, no strong gradients), the accuracy of those measures shows that 
current measurement techniques, applied to high-quality data, are reliable enough to allow 
a good determination of large-scale lensing statistics. 

3 Varying the S/N ratio on the shear map 

In the weak lensing approximation one usually estimates the shear inside a window by the 
mean of the galaxies' ellipticities inside the window, i.e. 
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So one considers the variance of the estimator 
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Figure 1: Variance and skcwncss of the convergence field measured from tire simulated survey (dash-dotted lines) compared to the theoretical 
values (solid lines). The "error" bars are the minimal errors calculated by van Wacrbeke ct al (1999). 



where n is the mean number of objects per window and O^^n means averaging on the 
elhpticity distribution and the spatial distribution of tracers. It is then clear that the signal- 
to-noise ratio of the shear depends on the local density of tracers, leading to local noise 
variation. To investigate this effect, we performed a simulation of a 25 deg^ survey, again 
with 2-dimensional shear maps; the positions of tracers were drawn either from a random 
(Poisson) distribution, or as the result of a random-walk process, of which third and fourth 
moment are independant of scale, while the angular autocorrelation function of the galaxies 
has an amplitude roughly one order of magnitude higher than that of faint blue galaxies (~ 
0.15 at 1', see e.g. Brainerd & Small (1998)). We used 60 realizations of the dark matter 
maps, and for each of them we have distributed intrinsically elliptical tracers {<7\e\ ~ 0-2), 
distributed in the two differents manners, and we studied 3 different densities (10, 30 and 
200 objects/arcmin^). 



The results on the variance for n — 30 galaxies/arcmin^ are show in fig. 2.; an excess of power 
at small scales is visible in the case of clustering of galaxies. Note however that in this set 
of simulations galaxy positions are uncorrelated with dark matter distribution (Bernardeau 
1998). Although this is not presented on a plot, we checked that the effect on the skewness 
is only due to the alteration of the variance. To correct for the effect, one has to realize the 
simple fact that the noise correction must not be done by means of eq. (|^) , but instead that 
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where < -i- > is filtered at the scale of interest. For strongly correlated, low-density samples, 
the < > statistic can significantly differ from 4. We can see on fig. 2. that this suffices to 
solve the problem. Such a correction is not as accurate in the n = 10 galaxies/ arcmin^ case, 
but this should not constitute a problem in an observational context as the current surveys 
have galaxies samples with larger density and correlation length. 
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FigUr6 2: The variance of the convergence held for a 25 sq degrees survey. Top panel : the tracers density is 30 objects/arcmin'^ and 

,,2 

a the ellipticity noise correction made by . Tire solid line is the theoretical value, the dotted line is the value obtained for a poisson 
distribution of tracers, the dashed-line is obtained when the tracers positions arc drawn from a random-walk. Bottom panel : the dashed line 
is the variance of the random-walk case, with noise corrected as cr^ X < > . The dash-dotted line is the same case with a tracer density of 
10 objects/arcmin'^ . For clarity only error bars of the 'problematic' cases are shown. 



4 Conclusion 

In the framework of future wide field surveys (e.g. the MEGACAM lensing survey), it it 
essential to verify the stability of the statistics of the density field through the measurement 
process. We checked this stability in two ways: first we performed a detailed simulation of 
a large (100 deg^), compact survey and analyzed the obtained images, and we succeeded to 
extract the cosmological signal. Secondly, we investigated the effect of rapidly varying S /N 
ratio on the convergence reconstruction via realizations of a random- walk process; the effect 
introduces spurious variance at small scales but can easily be corrected for in most relevant 
cases. This work therefore enhances our confidence on the capability of future wide surveys 
to accurately measure cosmic shear statistics. 
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